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Abstract. In order to extend the concept of geometric phase to discrete dynamics, we 
investigate its presence in discrete adiabatically-rotated rotators. For the rotated sine circle 
map, we demonstrate an analytical relationship between the geometric phase and the rotation 
number of the system. For the discrete version of the rotated rotator considered by Berry, the 
rotated standard map, we further explore this connection as well as the role of the geometric 
phase at the onset of chaos. Further into the chaotic regime, we show that the geometric phase 
is also related to the diffusive behavior of the dynamical variables and the Lyapunov exponent. 
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1. Introduction 

The geometric phase 0] is an example of anholonomy: the failure of certain variables to 
return to their original values after a closed circuit in the parameters. Physical expressions of 
such anholonomies appear in the rotation of the plane of oscillation of a Foucault pendulum 
ED, how swimming is performed by microorganisms at low Reynolds numbers [3J, and how 
a falling cat can manage to reorientate itself in mid air in order to land on its feet [4 J . The 
geometric phase was originally encountered — as Berry's phase — in quantum mechanics 
||5l . From there, it was generalized to classical integrable systems as Hannay's angle (6]|. 
Later it was extended to nonintegrable perturbations of Hamiltonian systems B71 [H, and 
thence to dissipative systems E QI1, all these instances within the context of continuous- 
time dynamics. In the same context, rotated rotators have been natural models in which to 
study this phenomenon because they provide an easy way to control the adiabatic nature of 
the cyclic variation of the parameter. With this in mind, Berry and Morgan, for instance, 
investigated the geometric phase of a continuous -time Hamiltonian rotated rotator [fTT||. 

In spite of the extensive research that has taken place in the last few decades on geometric 
phases in a large class of applications, neither the geometric phase nor any of its cognates have 
been considered hitherto in discrete dynamical systems. Our purpose in this Letter is to make 
good this deficit and to show that the geometric phase is linked to important aspects of the 
dynamics of maps. In order to do so, we adopt the paradigm of the rotated rotator but follow 
an inverse sequence to the historical development. In the first place we deal with the sine 
circle map that may be thought of as the discretization of a kicked rotator of the type Berry 
and Morgan studied with the addition of strong dissipation. We consider the results of discrete 
adiabatically-evolving parameter loops in such a prototypical discrete-time dynamical system. 
The geometric phase in the rotated circle map, it turns out, is intimately related to the behavior 
of the rotation number of the map as a function of the bare frequency parameter. Turning to 
the Hamiltonian side, we study the rotated standard map, in which we discover surprising 
relationships between the geometric phase, not only with the rotation number as in the former 
case, but also with the Lyapunov exponent and the diffusive behavior of both action and phase 
variables. The reason for our following this inverted developmental sequence is that ID circle 
maps are simpler as regards the transition from integrability to chaos, in comparison with the 
much richer behavior of the 2D Hamiltonian case. 

2. The rotated circle map: 

Before introducing our rotated version of the circle map, let us first recall a few necessary 
definitions and results on the original non-rotated one. The circle map, usually written as 

0n+i = faf{0o) = 9 n + n-^ sin27T0 n mod 1, (1) 

where 9 n = /^(0 O ) represents the nth iterate of 6q, which qualitatively describes the 
dynamics of two interacting nonlinear oscillators, is a one-dimensional discrete mapping that 
describes how a rotator of natural frequency f2 behaves when forced at frequency one through 
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Figure 1. (a) Devil's quarry plot of rotation number p in the sine circle map; a section through 
the quarry with k constant is a devil's staircase, (b) Geometric phase <fi g plotted against 57 for 
the rotated critical (k — 1) sine circle map of Eq. Q. (c) In a colormap, same as in (b) but as 
a function of both parameters f2 and k. 



a coupling of strength k. When k = the rotator runs uncoupled at frequency f2, but when 
k > it can lock into a periodic orbit: a resonance with some rational ratio p/q to the driving 
frequency. To measure the frequency of the rotator, i.e., the average rotation per iteration of 
the map, it is useful to define the rotation number 

On — 9o 



P 



lim 

AT-s-oo 



N 
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If we plot rotation number p against Vt and k — a few hundred iterations after discarding 
an initial transient are sufficient to give an accurate value for p — we obtain the devil's 
quarry illustrated in Fig. [jja)- Periodic orbits with different rational rotation numbers show 
up as so-called Arnold tongues: flat steps in Fig. [TJ a )- When k < 1 the map is termed 
subcritical, and intervals on which the rotation number is constant and rational, where there 
is a periodic orbit of a particular period, punctuate intervals of increasing rotation number, 
whereas in supercritical circle maps (k > 1) the periodic orbits overlap. Chaos is found 
in the supercritical circle map as iterates wander between the overlapping resonances. In a 
critical circle map at A; = 1, at every value of f2 there is a periodic orbit, and the rotation 
number increases in a staircase fashion with steps at each rational rotation number and risers 
in between. The devil's quarry becomes the devil's staircase when we look at a section with 
k constant through the quarry. The ordering of periodic orbits in the devil's staircase has 
been understood in terms of Farey sequences and Stern-Brocot trees ffT2l IT3l IT4l fT3Tl . and the 
transition to chaos in the circle map is well understood. 

Now let us consider the rotated circle map; probably the simplest discrete time system 
with a discretely and adiabatically varying parameter. To this end we introduce a discrete 
slowly varying parameter X n : 

k 

9 n+ i = 8 n + - — sm2ir(6 n + X n ) modi, 

Z7T 



X„ 



= X n + AX = X n ± 1/N, (3) 

where n = 1, 2, ... N, with N — » oo for adiabaticity. Is there a geometric contribution to the 
phase after such an excursion? Let us perform the change of variable 6' n = 6 n + X n , under 
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6' n+1 =e' n + (Q + AX) - — sin 2t< mod 1 



which the map can be written as 

k 

= 9' n + Q! - — sin 2n9' n mod 1, (4) 

where Q' = VI ± 1/N. So it is seen that the effect of the parameter loop is just a shift in the 
value of f2. 

In general, if one takes a system through a parameter loop, one obtains as a result three 
phases: a dynamic phase, a nonadiabatic phase, and a geometric phase. If one then traverses 
the same loop in the opposite direction, the dynamic phase accumulates as before, while the 
geometric phase is reversed in sign. There is, of course, still the nonadiabatic phase too; to 
get rid of this one must travel slowly around the loop. Thence the geometric phase may be 
obtained as 

<l>9 = j im ^ + n ^" > (5) 

where and 0_ are, respectively, the total angles 9 accumulated by traveling around the loop 
in the positive and negative directions. In terms of the primed variables, we can also define 

<f) g = hm ■ , (6) 

the obvious relation (j) g — cj) g ' — 1. Let us evaluate this limit, first for the simple case k — 0. 
Then /£ = 6 + NQ, so f£_ 1/Nfi = 9 + Ntt - 1 and f£ + y N>0 = 9 + NQ + 1, hence 
4> g ' = 1 and (f) g = 0. This limiting case is conceptually equivalent to that of a Foucault 
pendulum located at the Earth's equator where the plane of oscillations remains fixed as the 
Earth rotates. 

More interesting is what happens when k ^ 0. From the definition of the rotation number 

9 = i im — at = i im PnA e o)> (7) 

where we are defining Pn k (9 ) = (fnk{@o) ~ @o)/N, we have that 

, , _ fn+i/N,k(®o) ~ fn-i/N,k(Qo) 

Pn+i/N,k(® ^ ~ Pa-i/N,k(®o) /a , 

= lim ; . (8) 

N^oc 2/N 

We may recognize this as a peculiar type of derivative of the rotation number with respect to 
f2 in which the limit iV — > oo of the definition of p is taken simultaneously with the limit 
Af2 — > 0, maintaining NAQ = 1. We may denote such a limit as Bp/dVt. When both limits 
commute, this derivative coincides with the usual one and 

so that we can say, in this sense, that the geometric phase for the rotated circle map is the 
derivative of the rotation number. This is obviously true within the tongues where the rotation 
number is constant and both derivatives zero (Fig.[ljb and c)) and consequently <f> g ' = and 
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(f) g = — 1. One might think of these cases as equivalents of a Foucault pendulum located 
at one of the Earth's poles where the plane of oscillations rotates one turn a day in the 
opposite direction (minus sign) to the Earth's rotation. Outside the tongues both derivatives 
have singularities arising from the peculiar dependence of p on £7; nevertheless, apart from 
convergence (or lack thereof) details, 0' clearly has the qualitative aspect of a derivative of 
the devil's staircase, as illustrated in Fig. |TJb). Fig. [ijc), of which Fig. [TJb) is the k = 1 
slice, shows also two more things: for k < 1 where quasi-periodic orbits are still measure- 
wise abundant the phase outside the tongues is rather regular; this is a consequence of the 
well known "trivial" scaling relations of the subcritical mode-locking structure [?]. In the 
supercritical region k > 1 the limit is not well-behaved and fluctuations in the geometric 
phase appear that correspond to the chaotic regions already noted in Fig.[TJa). 

3. Berry's rotated rotator and the standard map: 

Berry [fTTTl considered the geometric phase of the rotated rotator q = V sin(g — X(t)), that is 
p = V sm(q-X(t)), 

q = P, (10) 

which is a rotator q = Vising being rotated by X(t). In the same way that discretizing 
the rotator with mixed forwards and backwards Euler methods gives us the area-preserving 
Chirikov-Taylor standard map, 

tv t o\ ■ J 7 "+! = In + ( fc / 27T ) M^0 n ) n i x 

I — °n ~r J n+1 

we can consider the rotated standard map that comes from the rotated rotator 

4+i = I n + (k/2ir) sin(27r(0 n + X n )), 

9n+l = ®n + 4+1 7 X n+ i — X n ± 1/N, (12) 

where X n is again a discretely and adiabatically varying parameter that is moved around a 
closed loop from to 1. 

In the rotating frame of reference, where 

n = n ± X n , (13) 

the rotated standard map reads: 

4+i = 4 + (k/2n) sin[(27T0 ± X n )}, (14) 
Ki+i -F X n+ \ — 9 n =p X n + 4+i, 

or 

9' n+1 =9' n ± (X n+1 - X n ) + 4+i = 9' n ±A + 4+1- (15) 
If we define I' n = I n ± A we obtain the original standard map for the new primed variables: 



4+i W 4 + (V27T) Bia[(20] 

9'n+l J I Q'n + 4+1 



(16) 
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Figure 2. The logarithm of the geometric phase for the rotated standard map in a rainbow 
color code as a function of the initial conditions (Iq,8o) for k = 1. (b) A transversal cut 
through (a) at 8q = 0.5 demonstrates that the phase is zero in the islands. 



By the definition of the geometric phase, 



lim 



and 



(17) 
(18) 
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[T"] g (Jo + 1/N, 9 ) - [T N ] e (Jo - 1/N, 9 C 
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Given that the rotation number p can be expressed as 

\T N ]Ji ,e )-e 

p Jo, O = I™ Piv (Jo, do) = Mm L Jel ? (19) 

we derive the following relation for the geometric phase: 

r Mio + 1/N, gg) - p jV (/ - 1/iV, gg) 

G>- = lim ; ; r 

= (20) 

Wo 

Of course, this limit does not necessarily exist, nor can one expect it to be finite, so 
we should take these equations as formal relations. However, as Fig. [2] shows, within the 
regular KAM islands where the rotation number remains constant, this relation is easily 
verified and the geometric phase is null. The situation is much less intuitive within the 
chaotic regions of the phase space; Fig. [2] again. There, the rotation number is not well 
defined and its dependence on the initial action has not been studied. Naively one could 
estimate the difference between two neighboring trajectories in the numerator of the r.h.s of 
Eq. ( 18 1 as proportional to 2 exp(AiV) /N, which clearly diverges exponentially. On the basis 
of this assumption one could attempt to establish a connection between this divergence of 
the geometric phase and the Lyapunov exponent A of the trajectory as In <p' g m NX + In N. 
However, the experimentally computed phase does indeed diverge with N, but at a much 
slower pace than exponentially. Instead, we find power-law behavior for (fi g with an exponent 
of 3/2, which is independent of the value of the nonlinearity parameter k (Fig. [3^). In other 
words, the phase difference of neighboring trajectories varying by an amount 1/N in their 
initial actions grows superdiffusively at the same rate as that of Richardson's law, after iV 
iterations. The failure of the exponential assumption is because the Lyapunov exponent only 
holds if the double limit of the initial perturbation going to zero and the time to infinity is 
taken in this order. The phase, however, is defined with a similar double limit but taken along 
the particular direction in which the product of time times the initial perturbation is kept equal 
to one. On the other hand, the 3/2 power law can be explained from the well-known diffusive 
behavior of the action due to the chaotic dynamics of the standard map. From the equation 
for the phase we may realize that the phase difference A9 between two trajectories evolves as 

(Ag) n+1 - (Ag) n = (A/) n+1 . (21) 

It follows, for large n, that 

d ^ 6 ^ ~ IA/L oc n 1 ' 2 => |Ag| oc n 3 / 2 . (22) 
dn 

But there is a further surprising result in the behavior of the phase illustrated in Fig. |3j5. By 
extrapolating the power law up to In iV = we obtain the prefactor of the hypothesized 3/2 
power law; this prefactor plotted as a function of the nonlinearity amplitude k behaves very 




Figure 3. The geometric phase for the rotated standard map diverges superdiffusively for 
large k with power-law exponent (3 ps 3/2. (b) The y-intercept log(^o) (filled circles) follows 
the logarithmic fc-dependence of the Lyapunov exponent, A (open circles) for large k. 



Geometric phases in discrete dynamical systems 



9 



similarly to the Lyapunov exponent plotted as a function of the same parameter, up to an 
additive constant ps 2. The extension, origin and consequences of this connection are under 
current investigation. 

In summary, we have shown that the study of adiabatic modulations in maps and, in 
particular, the extension of the concept of geometric phase to this type of systems opens up a 
new avenue to understand their dynamical features. 
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